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We show that for ultracold magnetic lanthanide atoms chaotic scattering emerges due to a com¬ 
bination of anisotropic interaction potentials and Zeeman coupling under an external magnetic 
held. This scattering is studied in a collaborative experimental and theoretical effort for both 
dysprosium and erbium. We present extensive atom-loss measurements of their dense magnetic 
Feshbach resonance spectra, analyze their statistical properties, and compare to predictions from a 
random-matrix-theory inspired model. Furthermore, theoretical coupled-channels simulations of the 
anisotropic molecular Hamiltonian at zero magnetic held show that weakly-bound, near threshold 
diatomic levels form overlapping, uncoupled chaotic series that when combined are randomly dis¬ 
tributed. The Zeeman interaction shifts and couples these levels, leading to a Feshbach spectrum of 
zero-energy bound states with nearest-neighbor spacings that changes from randomly to chaotically 
distributed for increasing magnetic held. Finally, we show that the extreme temperature sensitivity 
of a small, but sizeable fraction of the resonances in the Dy and Er atom-loss spectra is due to 
resonant non-zero partial-wave collisions. Our threshold analysis for these resonances indicates a 
large collision-energy dependence of the three-body recombination rate. 

PACS numbers: 03.65.Nk, 34.50.-s, 05.45.Mt 


I. INTRODUCTION 

Anisotropic interactions are a central and modern tool 
for engineering quantum few- and many-body processes 
[T]. A prominent example of such an interaction is the 
long-range dipole-dipole interaction (DDI) acting for in¬ 
stance between polar molecules [2], Rydberg atoms [3], 
or magnetic atoms [4]. Over the years, fascinating quan¬ 
tum effects of the anisotropy have been observed, such 
as the d-wave collapse of a dipolar Bose-Einstein conden¬ 
sate [5] , the deformation of the Fermi surface [6] , and the 
control of stereodynamics in dipolar collisions [7]. More¬ 
over, the DDI is expected to give rise to a plethora of 
few- and many-body phenomena, which still await obser¬ 
vation, such as universal few-body physics BM , rotonic 
features [Mill], two-dimensional stable solitons [T2] , and 
the supersolid phase [T5] . 

Recently, atomic species in the lanthanide family be¬ 
came available to the field of ultracold quantum gases. 
The interaction between magnetic lanthanide atoms. 
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such as Er muni and Dy mm, is highly anisotropic. 
This is not only due to the long-range DDI, originat¬ 
ing from their large magnetic moment, but also to the 
shorter-ranged van-der-Waals interaction [I8j . which ex¬ 
hibits anisotropic contributions arising from the large or¬ 
bital angular momentum of their valence electrons. 

For magnetic lanthanides, which also include the suc¬ 
cessfully laser-cooled elements Ho m and Tm [20] . 
the orbital anisotropy is a consequence of a partially 
filled submerged 4f electron shell that underlies a closed 
outer 6s shell. This leads to an electronic ground 
state with a total atomic angular momentum / with 
j1. Consequently, in collisions between such atoms 
there exist (j -|- 1)^ non-degenerate (gerade) molecular 
Born-Oppenheimer (BO) potentials and a correspond¬ 
ingly large manifold of collision channels with associated 
molecular bound states. This is in sharp contrast to 
the one or two BO potentials encountered in alkaline- 
earth and alkali-metal atom collisions. In addition, the 
anisotropy or orientation dependence of the BO poten¬ 
tials strongly mixes collision channels with large relative 
orbital angular momentum (. between the atoms even for 
our ultracold collisions with a f = 0, s-wave initial chan¬ 
nel. The complexity of the molecular forces are reflected 
in a dense spectrum of Fano-Feshbach resonances as a 
function of magnetic field B, as recently observed in Er 
[Ill El] and Dy [22]. In Er a statistical analysis of the 
spacings between resonances has shown correlations that 
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revealed chaotic scattering. The data set of the initial 
Dy experiments was too small to extract statistically- 
significant correlations. 

Chaotic behavior is manifest in a variety of complex 
systems ranging from atomic to nuclear and solid-state 
physics. In atomic physics chaos was originally studied 
with Rydberg states of H and He in a magnetic field 
[23]. Later on a variety of more complex atoms and 
ions in highly-excited states have shown signatures of 
chaotic spectral distributions |21|. The origin of chaos 
in these systems was traced back to a strong mixing of 
many-electron excited states by the Coulomb interaction 
[25| . A chaotic level distribution is also common in a 
variety of solid state systems ranging from those with 
strong many-body interactions to the motion of parti¬ 
cles in irregular potentials [MIEI]. Experiments in nu¬ 
clear physics US] |5S| have also produced substantial ev¬ 
idence for chaotic neutron resonance spectrum fluctua¬ 
tions, which agree with predictions of random matrix the¬ 
ory (RMT). Similar agreement was found from numeri¬ 
cal simulations based on nuclear shell models (SOI ED- 
Moreover, Refs. [311|3S] suggested that chaos is a generic 
property of nuclei with multiple degrees of freedom (i.e. 
multiple active shells), which become completely mixed. 

This article describes a joint effort to understand ul¬ 
tracold scattering and Fano-Feshbach spectra of strongly 
magnetic Er and Dy atoms. In particular, we report 
on the measurement and statistical analysis of Fano- 
Feshbach spectra for Dy and Er between B = 0 G to 
70 G at gas temperatures T below and around I fiK. Here 
1 G=0.1 mT. We observe that both elements have similar 
chaotic scattering. We present a random-matrix-theory 
(RMT) inspired model to gain insight into their statis¬ 
tical properties as well as theoretical evidence based on 
coupled-channels calculations with a microscopic Hamil¬ 
tonian that chaotic scattering requires both strong molec¬ 
ular anisotropy and Zeeman mixing to fully develop. 
Limitations of the RMT are also discussed. Finally, we 
present experimental data and a comparison to a reso¬ 
nant trimer model to show that our increase in resonance 
density with temperature is a consequence of the strong 
collision-energy dependence of transitions from entrance 
d-wave channels of three free atoms to resonant trimer 
states. 


II. EXPERIMENT 

A. Measurement 

The experimental study of Fano-Feshbach resonances 
in Er and Dy is based on high-resolution trap-loss spec¬ 
troscopy on spin-polarized thermal samples. Ultracold 
bosonic ^®'^Dy samples are created by direct loading from 
a narrow-line magneto-optical trap (MOT), operating on 
the 626 nm cycling transition, into a single-beam optical 
dipole trap (ODT) [Slj. By moving the last focusing lens 
of the ODT, the atoms are transported from the MOT 


chamber to the science cell. The ODT is created with a 
100 W fiber laser at a wavelength of 1070 nm. We achieve 
a transport efficiency close to unity. This fiber laser, how¬ 
ever, causes atom loss due to its longitudinal multimode 
structure |35j . Therefore, we transfer the atoms into a 
second single-beam ODT, created by a 55 W solid-state 
laser at a wavelength of 1064 nm. Finally, forced evap¬ 
orative cooling in a crossed ODT leads to a sample of 
10® atoms in the energetically-lowest Zeeman sublevel 
ruj = —8 at T = 600 nK. 

High-resolution trap-loss spectroscopy is performed on 
a spin-polarized bosonic ^®®Er sample at T = 1400 nK 
and compares this spectrum with that obtained at a four 
times lower temperature measured in previous work both 
for ®®®Er as well as for fermionic ^®^Er [21]. The ex¬ 
perimental procedures for creating bosonic and fermionic 
samples are described in Ref. m and ng, respectively. 
Bosons (fermions) are prepared in the lowest Zeeman 
sublevel rrij = —6 (m/ = —19/2). Erbium samples are 
trapped in a crossed ODT and contain about 10® atoms. 


B. Feshbach spectroscopy 

Feshbach spectroscopy is performed in a similar man¬ 
ner for the two species. The magnetic field is ramped 
up over a few milliseconds to a magnetic field value H, 
where the atoms are held in the ODT for 500 ms for Dy, 
400 ms for ^®®Er, and 100 ms for During this time 

inelastic three-body recombination causes atom loss from 
the ODT. At resonance, the recombination process is 
enhanced because of the coupling between the atomic- 
threshold state and a molecular state leading to a reso¬ 
nant increase of the atom loss. We identify the field loca¬ 
tions of maximum loss as the positions of Fano-Feshbach 
resonances [36j . The atom number is probed by standard 
time-of-flight absorption imaging at low magnetic field. 
We record atom-loss features for magnetic field values 
between OG and 70 G in steps of a few mG. Figure a) 
shows the normalized loss spectrum for the ^®^Dy iso¬ 
tope, where we identify 309 resonances. For ^®®Er at 
T = 1.4/rK there are 238 resonances. The Fano-Feshbach 
scan of fermionic ^®^Er is carried out from 0 G to 4.4 G 
and yields 115 resonances. 

The understanding of the richness of the scattering in 
Er and Dy requires the development of sophisticated mi¬ 
croscopic coupled-channels scattering models. We defer 
such analysis until later in this paper and first analyze 
our data following the statistical approach based on the 
random-matrix theory (RMT) advocated by Ref. [2T]. In 
particular, we study the correlations between resonance 
locations via the nearest-neighbor spacing (NNS) distri¬ 
bution and setup a RMT-like model, which accounts for 
the structure of our H-dependent microscopic Hamilto¬ 
nian to get intuition about these NNS. In our description 
of the coupled-channels calculations limitations of such a 
RMT-like model are discussed. 
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FIG. 1. (color online) (a) Trap-loss spectroscopy mapping of the Fano-Feshbach spectrum of as a function of magnetic field 

between B =0 G to 70 G with a data point every 14.5 mG and temperature T = 600 nK. Each data point is an average of three 
measurements, (b) Staircase function for the number of resonances as a function of B for ^®‘^Dy, ^®®Er at two temperatures, and 
fermionic ^®^Er. Dashed lines are linear fits forced to pass through the origin. Their slopes give a mean density of resonances 
of p = 4.3 G"^ for 2.7 G"^ for i®®Er at T = 350 nK, 3.4 G“^ for Er ^ T = 1.4/rK, and 25.6 G"^ for ^^'^Er. (c) Number 

variance of the experimental data as a function of scaled B-field interval N = AB]}. The experimental data for '^^^Dy (blue 
line) and ^®®Er at T = 350 nK (orange line) lies between the variances for an uncorrelated Poisson distribution (dashed line) 
and the correlated Wigner-Dyson distribution (dot-dashed line). 


C. Statistical analyses 

Our statistical analysis starts with the construction of 
the staircase function, which is a step-like function that 
counts the number of resonances below magnetic-field 
value B m- Figure [3b) shows the staircase function 
for Dy and Er. For both species the function is well-fit 
by a linear curve forced to pass through the origin. Its 
slope p corresponds to the density of resonances. Devi¬ 
ations below and above the fit occur for small and large 
B, respectively. The fitted resonance densities are given 
in the caption. Remarkably, the density of resonances of 
^®®Er at T = 1.4 pK is 25 % higher than the one observed 
at 350 nK. The discussion of the origin of this sensitivity 
is postponed until Section jV] The density p for bosonic 
Dy is 50% larger than for bosonic Er. This is caused by 
the larger f of Dy and, thus, its larger number of allowed 
collision channels. The much larger density p of 25.6 
for the fermionic ^®^Er is due to its additional hyperfine 
structure. 

Fluctuations in the number of resonances within a 
magnetic held interval AB is a second measure of the 
statistical properties of the spacings between resonances. 
Formally, it is dehned as the dimensionless number vari¬ 


ance = iV2 - where N = N,/M, = 

and Ni is the number of resonances in 
held interval [iAB, {i + l)Ai?] with i = 0,..., M — 1, 
such that MAB = B^^ax and Bmax = 70 G for both 
species. Consequently, N = ABp. For shot noise or a 
Poissonian distribution we expect = N. Figure j^c) 
compares for our Dy and Er data as a function of N. 
The huctuations for both species monotonically increase 
with AB but are substantially less than the shot noise 
limit. While this behavior was previously demonstrated 
for Er [21], the present results provide the hrst evidence 
of correlation in Dy and indicate similarity between the 
species. 

These correlations between resonance locations are fur¬ 
ther studied using the nearest-neighbor spacings distri¬ 
bution P{s) where s = SBp and 6B is the held spacing 
between two adjacent resonances in the spectrum. Fig- 
uresj^a) and (b) show the computed NNS distribution of 
our experimental data, derived from the number of NNS, 
Si, with spacings s between iSs and (i + l)(js, where 
j = 0,1,... and Ss « 0.3. The NNS distributions have 
clear deviations from both the Poisson Pp(s) = exp(—s) 
and Wigner-Dyson Pwd(s) = (7r/2)sexp[—(7r/4)s^] dis¬ 
tribution, two well-known distributions within RMT m- 
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FIG. 2. (color online) (a) Nearest-neighbor spacing distribution, P{s), determined from all-observed Fano-Feshbach resonances 
for (blue markers). Dashed and dot-dashed curves are the Poisson and Wigner-Dyson distribution, respectively. The 

solid line is a Brody distribution with rj — 0.45(7) ht to the experimental data, (b) Distributions P{s) for ^®®Er at T = 350 nK 
and T = 1.4pK (orange and red markers, respectively), and ^®^Er at T = OATf (green), where Tp is the Fermi temperature 
of the gas. The solid line is a Brody distribution with rj = 0.68(9) fit to the ^®®Er data at T = 350 nK. Panels (c) and (d) show 
the magnetic-field-resolved Brody parameter ri(B) as a function magnetic field for ^®"^Dy and ^®®Er, respectively. The Brody 
parameters for the Poisson and Wigner-Dyson distribution are 0 and 1, respectively. Gray markers and lines in all four panels 
are results from our coupled-channels calculations. The Icr error bars in (a) and (b) correspond to Poisson counting errors 
while shaded bands in (c) and (d) are Icr statistical uncertainties of the fits to the data. 


A Poisson distribution corresponds to a random distribu¬ 
tion of resonance locations, while a Wigner-Dyson dis¬ 
tribution corresponds to a situation, where neighbor¬ 
ing resonances “avoid” each other and Pwd(s) oc s for 
s —>■ 0. Deviations are also seen in Fig. [^c), where 
for both atomic species the variance T? does not agree 
with the corresponding predictions for these distribu¬ 
tions. The experimental NNS distributions in Figs, [^a) 
and (b) have also been fit to the Brody distribution 
Pb(s, rf) = 6(1 + ri)s'^ exp[—an empirical function 
that interpolates between Pp(s) and Pwd(s) for 77 = 0 
and 1 , respectively, and 6 is a normalization constant 
|38| . The values for 77 reported in the caption, indicate 
intermediate or mixed behavior of the data. 

We present the magnetic-field resolved Brody param¬ 
eter rj{B) in Figs, [^c) and (d) obtained from a fit to the 
NNS distribution of resonances located in moving inter¬ 
vals [B — AB/2,B + AB/2] with AB = 20 G. It has a 
non-negligible Icr uncertainty equally limited by the qual¬ 
ity of the fit and the number of Feshbach resonances in 
an interval or bin. The latter uncertainty is reflected in 
the bin-to-bin variation of ri{B). For Dy we observe that 
77 increase linearly with field for small S, which saturates 
at a value of « 0.5 for B > 30 G. For Er the Brody pa¬ 
rameter fluctuates around 0.5. Interestingly, the Er data 


at our two temperatures have a similar behavior, indi¬ 
cating that the larger density of resonances at higher T 
does not impact the degree of correlation between their 
spacings. 

III. RMT ENSEMBLE MODEL 

Random matrix theory is based on the powerful notion 
that the statistics of eigenvalues and eigenfunctions of a 
complex system can be studied by replacing the micro¬ 
scopic Hamiltonian by an ensemble of random Hamilto¬ 
nians. In this spirit, we construct a RMT-inspired model 
for weakly-bound molecular dimer states to test the dis¬ 
tribution of Eano-Eeshbach resonances. 

Our RMT model is based on the statistics of eigen¬ 
values of the N X N real, symmetric matrix F7 rmt = 
Hq + TJz, where matrices Hq and Hz represent the 
B = 0 Hamiltonian and the Zeeman interaction of the 
two atoms, respectively. Without loss of generality we 
can assume that Hz is a diagonal matrix with matrix 
elements given by rngfisB, where m is an integer be¬ 
tween —2j and 2 j, corresponding to the sum of the 
projection quantum numbers of the atomic angular mo¬ 
menta, g is the atomic Lande factor, and 773 is the Bohr 
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magneton. The Zeeman interaction does not depend on 
the rotational state of the molecule and, thus, entries 
in Hz correspond to states with a definite value for i 
and its projection. Hq is then the B = 0 Hamiltonian 
expressed in this basis. It is also convenient to define 
Hq = iJd + Hcpi, where diagonal matrix contains the 
diagonal matrix elements of Hq and H^pi is the matrix of 
all its off-diagonal elements. The eigenvalues of Hq can 
then be interpreted as the energies of ro-vibrational levels 
of the isotropic contribution of the molecular BO poten¬ 
tials, while i?cpi describes mixing due to the anisotropic 
contributions of these potentials. 

We generate members of our ensemble of I/rmt by 



FIG. 3. (color online) Molecular spectrum and NNS distri¬ 
butions of Fano-Feshbach resonances of ^®®Er calculated from 
our RMT model, (a) Example of a spectrum of molecular 
binding energies as a function of B for Vcpilh — 2 MHz, 
Tjd = 0. Here, h is Planck’s constant. Squares aX E/h = 0 
MHz indicate crossings of molecular levels with the thresh¬ 
old of two rrij = —6 atoms and correspond to the position 
of Feshbach resonances. For the sake of visibility we only 
show the spectrum between i? = 0 G to 10 G. (b) NNS dis¬ 
tributions of simulated Feshbach resonances for rjd = 0 and 
Vcpilh = 0 MHz (circles), 2MHz (squares), and 10MHz (tri¬ 
angles). The dashed and dashed-dotted lines are Poisson and 
Wigner-Dyson distributions, respectively. Solid lines are best- 
fit Brody distributions with rj = 0.03(1), rj = 0.41(5) and 
rj = 0.82(1) for Vcpijh = 0 MHz, 2 MHz, and 10 MHz, respec¬ 
tively. (c) Fitted Brody parameters of the NN energy-spacing 
distribution of the eigenvalues of Hrmt at H = 0 G as a 
function of Vcpi- Gircles, squares, down triangles, and up tri¬ 
angles correspond to rfd = 0, 0.25, 0.5, and 1, respectively, 
(d) Fitted Brody parameters of the NNS distribution of the 
Feshbach resonances as a function of Vcpi and for four r]d using 
the same marker code as in (c). In (b), (c), and (d) distribu¬ 
tions are obtained by averaging over 15 realizations of Hrmt, 
each of dimension 500 x 500 and using Feshbach resonances 
computed up to B = 85 G. 


choosing random matrix elements for Hz, Hq, and i?cpi 
based on specific distributions. The values of m in Hz are 
uniformly distributed integers between —2j and 2j. The 
matrix elements of Hq are chosen according to a Brody 
distribution with variable Brody parameter r]Q € [0,1] 
and with a mean energy spacing between bound states ea. 
Finally, matrix elements of H^pi are chosen as Gaussian- 
distributed real numbers with zero mean and standard 
deviation v^pi, thereby on average coupling all diago¬ 
nal elements equally. Notice that this construction devi¬ 
ates from that for a true Gaussian orthogonal ensemble, 
where all matrix elements of a symmetric Hamiltonian 
are Gaussian distributed [33] ■ 

We apply the RMT model to the case of ^®®Er. The 
relevant species-specific quantities are j = 6 , g = 1.16, 
and Cd is chosen to roughly reproduce the observed den¬ 
sity of Fano-Feshbach resonances of ^®®Er and is set to 
^d/h = 6.4 MHz, where h is Planck’s constant. Fig¬ 
ure l^a) shows an example of a molecular spectrum, the 
eigenvalues of FIrmt obtained with our RMT model as a 
function of B with ? 7 d = 0 and Vcpxjh = 2 MHz. We ob¬ 
serve that as B increases weakly-bound molecular states 
avoid each other multiple times before reaching the two- 
atom threshold creating a Feshbach resonance. When we 
turn off Hcp\ the levels cross. Similar H-field dependen¬ 
cies of the eigenvalues occur for r]Q > 0 . 

We investigate the effect of the parameters r'cpi and Vd 
on the NNS distribution of the Fano-Feshbach resonances 
as well as that of the B = 0 molecular levels. Figure 
l^b) shows the NNS distribution of Feshbach resonances, 
obtained by averaging over 15 realizations of FIrmt, for 
four values of v^px and rjQ = 0. For negligible v^px the 
distribution follows Fp(s) and approaches Fwd(s) when 
the anisotropic coupling strength z^cpi is large compared 
to Ed. In fact, we find that a larger Ed requires a larger 
Vcpi to develop correlations. 

Figuresj^c) and (d) show Brody parameters fit to NNS 
distributions as functions of r'cpi and rjQ. Panel (c) shows 
r] for the B = 0 molecular binding energies. For z/cpi = 0 
the Brody parameter is simply rjQ as expected from the 
distribution of the diagonal Hq, while for larger inter¬ 
action anisotropy Vcpx the parameter rj « 0.9, close to 
a Wigner-Dyson distribution, independent of tjq. Panel 
(d) shows r] extracted from the RMT Feshbach resonance 
locations as a function of r'cpi- It suggests that the cor¬ 
relation in the NNS of the resonances is caused by z/cpi, 
whereas it appears fairly independent on r]Q. More pre¬ 
cisely, the Brody parameter fitted to these distributions 
rapidly increases from ry « 0 to 77 « 0.8 for z/cpi ^ Cd 
and tends to one for larger z/cpi- We conclude from the 
RMT model that the correlations between the locations 
of the Fano-Feshbach resonance are essentially due to the 
avoided crossings between weakly-bound molecular states 
at finite B and only weakly dependent on the energy dis¬ 
tribution at H = 0. In fact, these correlations increase 
for increasing z^cpi. 
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IV. MICROSCOPIC COUPLED-CHANNELS 
MODEL 

A. Realistic setup 

A quantitative understanding of the origin of the 
chaotic resonance distribution requires coupled-channels 
and bound-state calculations with physically realistic 
angular-momentum couplings and interaction potentials. 
We do so here based on the time-reversal symmet¬ 
ric Hamiltonian for the relative motion of Dy and Er 
described in Refs. HU EO]. It contains the Zeeman 
Hamiltonian, the molecular vibration and rotation, and 
the molecular interactions with isotropic (orientation- 
independent) and anisotropic (orientation-dependent) 
contributions, Vi{R) and Va[R), respectively, where R 
describes the separation R and orientation of the atom 
pair R. The potential has eight tensor operators coupling 
the two atomic and relative orbital angular momenta, 
Ji, j 2 , and For B = 0 the total angular momentum 
J = ji + j 2 + i IS conserved. For B > 0 G only the pro¬ 
jection M of J along B is conserved. The zero-of-energy 
of the Hamiltonian is the energy of an atom pair in the 
absolute lowest Zeeman sublevel rrija = —ja- 

The potentials Vi{R) and Va,{R) contain short-ranged 
exchange, medium-ranged van der Waals as well as long- 
range magnetic dipole-dipole interactions. We use the 
isotropic van der Waals coefficient Cg = 1723F;iaQ and 
anisotropic coefficients spread over ACg = ITAEhaQ for 
Er [21j . For Dy we have improved the value of van 
der Waals coefficients of Ref. [40] by including addi¬ 
tional experimental and theoretical transition frequen¬ 
cies and oscillator strengths mHn and now use Cg = 
2003 and spread ACg = 188 Ehag. In particular, 
the anisotropic spread for Dy has significantly increased. 
Here, Ef^ = 4.360 x 10“^® J is the Hartree energy and 
oq = 0.05297 nm is the Bohr radius. 


B. Bound-state calculations 

In Ref. [2T] we performed initial coupled-channels cal¬ 
culations of the scattering between ultracold Er atoms 
and predicted that tens of partial waves £ should have 
been included as the strength of the anisotropic contribu¬ 
tion is large. We, however, were unable to reach numer¬ 
ical convergence with respect to the number of coupled 
equations. 

Here, we circumvent this limitation by performing 
multichannel bound-state calculations, in which we use 
5 = 0 eigenstates as a basis for those at 5 > 0 G. For 
5 = 0, where J is a good quantum number, at most 
49 and 81 Bose-symmetrized and parity-conserving chan¬ 
nels are coupled for Er and Dy, respectively. The 5 = 0 
coupled Schrodinger equations are discretized on the in¬ 
terval 5 G [0, 5max] assuming zero boundary conditions 
and solved as a matrix eigenvalue problem H5IH5] . For 


each J only eigenstates with energies between [5o,5i] 
surrounding the zero of energy are computed and stored. 
The bound states for 5 > 0 G are solutions of the matrix 
eigenvalue problem that includes all computed 5 = 0 so¬ 
lutions with \M\ < J < Jmax and their coupling due to 
the Zeeman interaction. Selection rules of the Zeeman in¬ 
teraction ensure that there only exists direct coupling be¬ 
tween J and J' zero-field eigenstates with J— J' = 0, ±1. 
For both species 5max = 1000og, E^/h = —3 GHz and 
Ei/h = 0.9 GHz ensuring that Feshbach resonance loca¬ 
tions below 70 G are converged. 

In this section on the microscopic calculations we fo¬ 
cus on analysing the spectra at our coldest temperatures, 
where the initial collision channel has s-wave {£ = 0) 
character. Hence, we only need to consider even-.^ chan¬ 
nels with total projection quantum number M = —12 
and —16 for ^®®Er and ^®^Dy, respectively and inclusion 
of zero-field solutions up to Jmax = 36 for Dy and 39 for 
Er is sufficient to reproduce the experimental resonance 
densities. In Sec. |V] we will discuss higher temperature 
collisions between Er atoms, where J-wave {£ = 2) en¬ 
trance channels must be considered and, hence, spectra 
at other M values (i.e. M between —14 and —10 for 
^®®Er) contribute. 

C. Interaction anisotropies 

We first look into the role of interaction anisotropies on 
the level distribution of the most-weakly bound molecu¬ 
lar energy levels at zero magnetic field. There are two 
dominant components to the anisotropy, the dispersion 
V/\Cs(.R) and magnetic dipole-dipole Umdd(-R) contribu¬ 
tion. To distinguish the contributions of these two terms, 
we define: 

K(.R) = ^ACe^ACaiR) + AmDdFmdd(.R) (1) 

with variable strength XaCb and Amdd- We systemati¬ 
cally increase the strengths XaCb and Amdd from zero, 
where we recover the full physical strength for Amdd = 
XaCb = 1- 

For completeness we note that the dominant tensor 
operator for the anisotropic dispersion contribution is 

Vac,{R) = 4 E ^ {3(^ • MR ■ 3^) + 

with strength < 0 found with the methodology dis¬ 
cussed in the subsection llV A1 Weaker contributions indi¬ 
cated by dots are included in our calculations. Moreover, 

14idd(.R) = ^^^3^ {^^(R ■ 3i){R ■ h) - Ji ■ Jij , 

where /ig is the magnetic constant. 

Figure l^b) and (d) show the most-weakly-bound 5 = 
0, J = 16 levels of ^®"^Dy 2 as a function of anisotropy 
strength for purely dipolar (Aac6 = varying Amdd) 
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FIG. 4. (Color online) Interaction anisotropy induced chaos of B = 0 near-threshold bound states, (b) Weakly-bound J = 16 
bound-state energies of ^®^Dy 2 as a function of the anisotropy scale Amdd with AaCb = 0- (s-) NNS distribution (red circles) 
for the J = 16 bound state data in panel (b) at Amdd = 1 and Xac^ = 0. The solid red line is a Brody distribution fit to the 
data and agrees well with a Poisson distribution, (d) Weakly-bound J = 16 bound-state energies of ^®^Dy 2 as a function of 
the anisotropy scale Aac^ with Amdd = 0. (c) NNS distribution (purple squares) for the J = 16 bound state data in panel 
(d) at Amdd = 0 and Aac6 = 1- The solid purple line is a Brody distribution fit to the data and is close to a Wigner-Dyson 
distribution, (e) Moving average of the Brody parameter 77 as a function of AaCb (purple squares) or Amdd (red circles) with 
bins AA = 0.2 obtained by fitting the NNS distribution for the J = 16 bound state data in panel (b) and (d) to Brody 
distributions, respectively. The horizontal lines at r; = 0 and 1 correspond to the Brody parameter for a Poisson and Wigner- 
Dyson distribution, respectively. The la error bars combine statistical and fitting uncertainties, (f) The individual-J (blue 
squares) and combined-J (red circles) NNS distributions P{s) at Amdd = Aac6 = 1 as a function of the normalized energy 
spacing s. The distributions are derived from 5 = 0 bound-state data for J = 16,..., 25. The grey shaded areas in panels (a), 
(b), and (f) indicate the Wigner-Dyson distribution. 


and dispersive (Amdd = 0 , varying AaCb) anisotropic 
interaction, respectively. For Amdd = AaCb = ^ 
binding energies are regularly structured with many near 
degeneracies. In fact, the corresponding states are ro- 
vibrational levels of the isotropic centrifugal potentials 
Vi (R) and labeled by £. In our 3 GHz energy window an s- 
wave channel has at most three bound states, while even 
£ > 0 channels with their centrifugal barriers have fewer 
pn |49]. For small Aac 6 Amdd the degeneracy is 
lifted and levels shift linearly. The linear dependence for 
increasing strength of the dipole-dipole is approximately 
valid up to the physical value of Amdd = 1- Hence, the 
dipole-dipole interaction does not lead to our chaotic level 
distributions. In fact. Fig. [^a) shows that at Amdd = 1 
and AaCs = 0 the NND distribution is Poissonian. 

On the other hand, for a relatively small anisotropic 
dispersion strength XaCb ~ O.I levels start to avoid each 
other. Starting from AaCg ~ t).5 most avoided crossings 
are noticeable on the 3 GHz scale of the figure. At the 
nominal Aac 6 “ where there are 56 levels with —3 
GHz < E/h < 0 GHz, a significant fraction of the levels 
have undergone multiple avoided crossings and can not 
be described by a single dominant partial wave. The level 


spacing is chaotic as confirmed by the NND distribution 
for Aac6 = 1 Amdd = 0 in Fig. [^c). We have 
computed the weakly-bound J = 16 levels for Amdd = 
AaC 6 = 1- Visually the level distribution is much the 
same as the one shown in Fig. |^d). Similar results have 
been obtained for ^®®Er 2 . 

Figure l^e) quantifies the intuition gained from panels 
(a)-(d) by showing the Brody parameter 77 of the B = 0 
J = 16 ^®‘^Dy 2 levels as a function of Xac^ or Amdd- The 
Brody parameter is obtained by fitting a Brody distribu¬ 
tion to the NNS distribution of the bound state data in 
Fig. I^b) and (d). For increasing dipole-dipole strength 
Amdd and no anisotropic dispersion (Aac 6 ~ 0) P^“ 

rameter is always zero indicating the prevalence of small 
level spacings. On the other hand, in the absence of the 
dipole-dipole interaction, increasing XaCs leads to an in¬ 
creasing 77 . It evolves from 77 = 0.2 for Aac 6 ~ 

77 = 0.7 for XaCb — indicating a depopulation of small 
energy spacings. Note that our systems does not reach a 
Wigner-Dyson distribution, which corresponds to 77 = 1. 

We compare in Fig. if) two NNS distributions of 
B = 0 weakly-bound states of ^®'^Dy 2 obtained for the full 
anisotropic interaction (Amdd = AaCs = !)■ Both distri- 
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FIG. 5. (Color online) Theoretical ^®^Dy 2 (panels (a) and (b)) and ^®®Er 2 (panels (c) and (d)) near-threshold bound states as 
a function of magnetic held. Calculations have been performed with the full nominal anisotropy. Panels (a) and (c) show the 
near-threshold region between B = 0 G to 10 G while panels (b) and (d) show the region between B = 50 G to 60 G. Red 
crosses indicate the location of Feshbach resonances. 


butions are based on \E/h\ < 3 GHz bound states com¬ 
puted for J = 16 up to 25. The first so-called individual- 
J distribution is constructed by averaging the NNS distri¬ 
bution of levels for individual Js assuming that individ¬ 
ual distributions are the same. The second, combined-J 
NNS distribution, is calculated from a sorted list of all 
J = 16,..., 25 levels. Data for J > 25 are not included 
as the number of bound states is too small for a reliable 
determination of the NNS distribution. 

The individual-J NNS distribution is non-Poissonian 
as levels with the same J repel each other. The 
combined-J distribution, however, follows a Poisson dis¬ 
tribution indicating that energies of bound states with 
different J are uncorrelated. In other words, even though 
the Hamiltonian, i.e. the set of coupling operators be¬ 
tween ji, J 2 and £, is the same, differences in the matrix 
elements and thus coupling strengths between channels 
lead to uncorrelated eigen energies. 


D. Atom scattering in a magnetic field 

The study of the B = 0 G multi-channel bound states 
has shown that interaction anisotropies mix channels 
with the same J, while states with different J are un¬ 
correlated. The Zeeman interaction mixes these molec¬ 
ular levels and leads to the Fano-Feshbach spectrum. 
Figures [^a) and (b) show example ^®"^Dy 2 M = —16 
bound- state spectra as a function of B on two binding- 
energy and field regions. Similarly, Figs, [^c) and (d) 
show M = —12 ^®®Er 2 bound states. In all cases the full 
nominal anisotropy (AaCb = ^ Amdd = 1 ) is used. 
For Dy and Er channels with J up to 36 and 39 are in¬ 


cluded, respectively. The figure shows that the Dy level 
density is higher than that for Er. This simply follows 
from the larger atomic angular momentum of Dy, leading 
to a larger number of channels with the same J — \M\. 
We also observe that for both species the level structure 
in the 0 G to 10 G, small field region is qualitatively dif¬ 
ferent from that in the larger field region. Eor small B 
the avoided crossings are substantially narrower than for 
larger B. Moreover, at small field the levels cluster, while 
at larger field they are more uniformly distributed. These 
changes are a consequence of the linearly-increasing Zee- 
man coupling between vibrational levels with different Js 
as a function of B. 

Eigure |^a) shows effective length as{B) as a func¬ 
tion of B. It diverges at every resonance location and is 
closely related to the scattering length of a zero-energy 
collision. Our calculations can not be directly used to de¬ 
fine the scattering length as we use a hard-wall potential 
for R > Bmax- This wall leads to a discrete set of states 
with positive energy and using the lowest of these Es{B) 
we can define the effective length as{B) shown in the fig¬ 
ure by solving for Bs(B) = ;i^7r2/[2^r(Bmax - ^^(B))^] 
with Hr = to/ 2 and atomic mass to [50j . 

It is of interest to briefly discuss the convergence prop¬ 
erties of our calculations. The data in Figs. I^a) and 
[^c) and (d) are based on computations with channels 
with J up to Jmax = 39. Figure [^b) shows the ^®®Er 2 
Feshbach resonance density p as a function of Jmax- The 
resonance density increase linearly from « 0.5 1/G at 
Jmax = 12 but then is seen to “saturate” for larger Jmax- 
At Jmax = 39 the experimental density is reproduced. In 
addition, Fig. 0 c) shows the field location of resonances 
between 50 G and 55 G as a function of Jmax- The res- 
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FIG. 6. (Color online) Theoretical near-threshold bound states and Feshbach resonances in a magnetic field, (a) Near-threshold 
bound states for M = —12 ^®®Er 2 for fields between B = 50 G and 60 G (bottom half of image). The top half of the image 
shows the effective scattering length, defined in the text, as a function of B. It is infinite at a resonance location where a 
bound state has zero energy. Calculations use the physical interaction anisotropies and channel states with J < Jmax = 39. 
(b) Theoretical Feshbach resonance density p as a function of Jmax (purple circles) computed from resonance locations between 
0 G and 70 G. The dashed horizontal line indicates the experimental density for ^®®Er 2 . (c) Convergence study of resonance 
locations (crosses) between 50 G and 55 G as a function of Jmax. Purple lines connect resonances when their location has 
converged. 


onance locations change significantly for Jmax < 22, but 
then rapidly converge. This implies strong mixing among 
bound states with those J. On the other hand, the lo¬ 
cation of resonances that appear for J > 22 are almost 
immediately converged indicating weak mixing to smaller 
J states. 


E. Comparison of experiment and 
coupled-channels model 

In Figs, [^a) and (b) we show the NNS distribution 
of converged Feshbach-resonance locations based on our 
multi-channel data between B = 0 G and 70 G for ^®^Dy 2 
with Jmax = 36 and ^®®Er 2 with Jmax = 39, respectively. 
For both species the distribution clearly deviates from 
a Poisson distribution, consistent with the experimental 
distributions that are also shown. The fitted experimen¬ 
tal and coupled-channel Brody parameters agree within 
their error bars. 

The anisotropy parameters AaCb Amdd in the 
coupled-channels calculations and the parameter zzcpi in 
the RMT play analogous roles in the Hamiltonian and in 
the emergence of chaotic level distributions, even though 


no explicit quantitative connection exists. This role is 
most manifest in the Brody parameters of the B = 0 G 
bound states and that of the Feshbach resonance spectra 
for the two models. For ^®®Er the corresponding Brody 
parameters from the coupled-channels calculations are 
« 0.01 and 0.68 at the physical Amdd = Aac6 = 
respectively. Within the RMT model the small rj value 
for the B = 0 G level distribution requires weak coupling 
i^cpi ^ Cd and ryd ~ 0. In contrast the Brody parameter 
for the Feshbach resonance spectrum requires v^-pi « ea 
and points at limitations of the current RMT model. 
Similar conclusions hold for bosonic Dy. Future advanced 
RMT models might circumvent these limitations by in¬ 
corporating overlapping, uncoupled chaotic series as is 
found from our B = 0 G coupled-channels calculations. 

We plot the B-held-resolved Brody parameter ry(B) of 
the theoretical coupled-channels data in Figs. I^c) and 
(d). A comparison with the experimental rj{B) shows 
excellent agreement for ^®^Dy, while the agreement for 
i68Er is less satisfactorily. A possible explanation for the 
discrepancies in ^®®Er is the larger bin-to-bin fluctuations 
as bins contain fewer resonances than for ^®‘^Dy. 

For ^®^Dy the theoretical held-resolved Brody param¬ 
eter in Fig. Sc) linearly increases from zero for small 
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B fields and saturates at r]{B) « 0.5 for fields larger 
than 35 G where the size or width of the avoided cross¬ 
ings between weakly-bound states is larger. For ^®®Er in 
Fig. I^d) we find a much more rapid increase of r]{B) at 
small fields. This is followed by a plateau at ri{B) « 0.5 
between B — 20 G and 50 G, after which ri{B) 0.9 
with an uncertainty of 0.2 close to a Wigner distribu¬ 
tion. The initial rise of r]{B) for both atomic species is a 
consequence of weakly-bound vibrational levels, uncou¬ 
pled and randomly distributed when B = 0 G, that start 
to repel each other as the Zeeman interaction increases 
in strength for increasing B. The plateau at r]{B) « 0.5 
and the sudden increase of ri{B) to one for ^®®Er have 
no simple explanation and are determined by the not- 
fully-explored complex interplay between the Zeeman 
and anisotropic inter-atomic interactions. It does, how¬ 
ever, indicate that Wigner’s assumptions on ensembles of 
Hamiltonians do not hold for fields below 50 G. 


V. TEMPERATURE DEPENDENCE OF THE 
RESONANCE DENSITY 

We now describe the origin of the strong temperature 
dependence of some of the resonances in our atom-loss 
spectra and thus explain the accompanying increase of 
the resonance density. Here, atom loss is solely due to 
three-body recombination, where three ultra-cold atoms 
collide to form a di-atomic molecule and an atom that 
both are lost from the atom trap. Figure l^a) shows 
atom-loss spectra for one such resonance for ^^Er at four 
temperatures below 2 ^K. Atom loss, indeed, is larger for 
larger temperatures, but we also observe a broadening of 
the H-field width and a shift of the maximum loss posi¬ 
tion to larger B fields. Resonances with a weak temper¬ 
ature dependence show none of these behaviors. 

We show with an intuitive resonant “trimer” model 
that a strongly temperature-dependent resonance is due 
to scattering processes with entrance d-wave channels 
even though the two-body d-wave centrifugal barrier, 
= 250 /J,K, is a hundred times larger than our 
highest temperature, where /cs is the Boltzmann con¬ 
stant. The difference in the power-law Wigner-threshold 
behavior of the recombination rate with collision energy 
for s- and d-wave entrance-channel collisions can explain 
our observations. 

Three-body recombination has been extensively stud¬ 
ied in the context of Efimov physics [5TH5^ . We follow 
Refs. [55H57| and start from a coupled-channels descrip¬ 
tion in the (mass scaled) hyperradius p, which describes 
the size of the three-atomic system, and basis functions 
in the five other hyperspherical coordinates that are p- 
dependent eigen states of the squared “grand-angular 
momentum operator”. Similar to the coupled-channels 
description for two atoms, there are entrance, open, and 
closed channels. The collision starts in one of the en¬ 
trance channels with atoms in the energetically-lowest 
Zeeman state and relative three-body kinetic energy 


the dimer plus atom are the open channels, and bound 
states in closed channels can lead to resonances. These 
closed channels dissociate to three free-atom states with 
at least one atom in a Zeeman level with higher internal 
energy. The bound states are resonant “trimer” states 
giving us our name for the model. It should, however, be 
realized that their origin lies in bound states of pairs of 
atoms and that the resonant state is better thought of as 
a pair bound state that hops from pair to pair. We define 
E 3 = ‘h?k 1 /{ 2 p 3 ) = P 3 V 3 I 2 with the three-body reduced 
mass p 3 = where fca and V 3 are the relative wave 

vector and velocity, respectively. 

The potentials in the entrance channels have long- 
range repulsive centrifugal potentials, governed by the 
asymptotic behavior of the grand-angular momentum op¬ 
erator and depend on the relative orbital angular mo¬ 
mentum M of the three atoms. In fact, the centrifugal 
potentials are ?i^(A -I- 3/2)(A -I- 5 / 2 )/( 2 /r 3 p^) with non¬ 
negative integer quantum number A [5S] . For M = 0 the 
least repulsive potential has A = 0, while that for Af = 2 
has A = 2. 

For an isolated trimer resonance in a closed channel 
coupled to both entrance and other open channels we can 
apply the resonance theories by Fano and Feshbach and 
derive that the recombination rate coefficient at collision 
energy E 3 and entrance channel with quantum number 
A is given by L 3 {E 3 ,B) = V 3 a{E 3 , B), where the cross 
section a{E 3 ,B) = {2N + l)1927r^|5'(E3,i3)p/A:| and 

iccp Ri |2 = _ r(£'3)Fbr _ 

{E3-p{B-B3)y + {V,M)m^ 

is a resonant expression for the square of a dimension¬ 
less S'-matrix element, where Bq is the trimer resonance 
location and p is the magnetic moment of the resonant 
trimer relative to that of the entrance channel. The defi¬ 
nition for \S{E 3 , B)\‘^ also contains the entrance-channel 
energy width r(£' 3 ) = AxE^^^ to the trimer resonance 
with a characteristic power-law energy dependence that 
reflects the threshold behavior of the scattering solutions 
in the centrifugal potentials. The energy width Fbr de¬ 
termines the decay or breakup rate of the resonance into 
the fast atom and dimer pair and is independent of E 3 . 
Finally, rtot(A 3 ) = ^{E^) + Fbr- For simplicity we have 
assumed that non-resonant, direct recombination from 
the entrance to open channels is weak. We also note 
that for Af = 0 and A = 0, ^ 3 ( 153 , B) approaches a finite 
constant for E 3 —>• 0 as expected. 

In our experiments we have thermal samples of Er and 
we require the thermally-averaged rate coefficient 

1 

Ls (T, B) = - E^dEL3{E,B) 

^ Jo 

and normalization Z = E^dE = 2{kT)^. In 

order to increase signal to noise we have allowed a signif¬ 
icant fraction of atoms to be lost (See Fig. [^c)), which 
assuming a homogeneous sample can be modeled by the 
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FIG. 7. (Color online) Line shapes for a strongly-temperature-dependent ^®®Er Feshbach resonance near B = 1.48 G. Panel 
(a) shows experimental data (markers with error bars) as a function of B of the remaining atom number divided by the atom 
number away from resonance measured 400 ms after initial preparation. Black, red, green, and blue markers correspond to 
data for temperatures T — 230 nK, 740 nK, 1400 nK, and 2000 nK, respectively. Dashed lines connecting the markers guide 
the eye. Solid lines are theoretical line shapes of the remaining atom number based on the d-wave (Af = 2) recombination rates 
shown in panel (c). Panels (b) and (c) show simulated three-body recombination rates for the same four temperatures assuming 
three-body entrance s- (Af = 0) and d-wave (A/” = 2) scattering, respectively. Curves are based on a thermally-averaged line 
shape discussed in the text. Recombination rates are scaled such that the largest value in each panel is one. For both panels 
fj, = S.lfiB and Pbr/A^s = 250 nK, while V{E)/kB = 0.2(if/i?ref)^ nK in panel (b) and V{E)/kB = 0.1(i?/i5ref)^ nK in panel 
(c), where E^ef/kB = 1000 nK. 


rate equation dn{t)/dt = —2>L‘i{T, B)v?{t) for atom den¬ 
sity n[t) [53] with solution 

iV(th, B) = , 

y/\ + QL^{T,B)nlt^ 

where B) is the remaining atom number after hold 

time th, Nq is the initial atom number, and no is the 
initial density. This non-linear time evolution adds addi¬ 
tional broadening to the lines. 

Figures Bb) and (c) show our model event rates 
L'j,{T^B) as a function of i? for Af = 0, A = 0 and 
AC = 2, A = 2 , respectively. Curves are for the same 
four temperatures as in Fig. j^a). A comparison of pan¬ 
els (b) and (c) shows a striking difference. The strongest 
features in panel (b) are for the smallest temperatures, 
while those in panel (c) are for the largest tempera¬ 
tures. This behavior naturally follows from an approxi¬ 
mation of the integrant in L 3 {T, B) under the conditions 
kT ^ Fbr r(if) [^. In this limit the Lorentzian is 
sharply peaked around = fi{B — Bq) for B > Bq and 
after some algebra it follows that L^iT^B) as a function 
of B has a maximum value proportional to {kT)^~^ lo¬ 
cated at B — Bq + {X + 2)kT j[i. Consequently, for A = 0 
and 2 the maximum loss rate decreases and increases with 
T, respectively. Even for less restrictive parameter values 
as used in Fig. [^this trend remains. 

Our experimental data has a temperature trend as in 
panel (c). In fact. Fig. I^a) compares our experimental 
loss data with model Af (ti,, B) for Af = 2, A = 2 using the 
same parameters as in Fig. j^c) and requiring a « 50% 


maximum atom loss as in the experiment. It is worth not¬ 
ing that, from our theoretical calculations, the magnetic 
field width of L^{T, B) is noticeably smaller than that for 
N{t\i,B) indicating that the finite hold time does indeed 
lead to broadening. The agreement of the experimental 
data and the prediction of our model for the losses is sat¬ 
isfactory for all four temperatures given the limitations 
and approximations within our modeling. We conclude 
that our strongly T-dependent resonances correspond to 
d-wave or more precisely Af = 2 entrance channel col¬ 
lisions. Note that we haven’t observed any resonances 
with temperature dependence similar to Fig. Bb) in our 
spectra. In the case of resonances with three-body s- 
wave entrance channel, which would correspond to such 
a dependence, we infer that the loss spectra is saturated. 
This will be subject for future investigations. 

As a corollary, this implies that for two colliding atoms 
as described in Sec. |IV| temperature-dependent reso¬ 
nances are due to collisions with entrance d-waves for 
which there are multiple allowed values of the total an¬ 
gular projection quantum number M. Here M = -14 to 
-10 for bosonic Er and M = -18 to -14 for bosonic Dy. 
Numerical computations, not presented here, show that 
their zero-field bound-states and thus resonance locations 
are again uncorrelated and random. 

VI. CONCLUSION 

In summary, we have experimentally and theoretically 
studied the resonant scattering of ultracold Er and Dy 
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atoms in a magnetic field. We have shown that chaotic 
scattering as witnessed by chaotic nearest-neighbor spac- 
ings between Feshbach resonance locations emerges due 
to the anisotropy in the molecular dispersion. 

Our study has also revealed several unique features of 
colliding magnetic lanthanides that have not been ob¬ 
served in any other ultracold atomic system. These lan¬ 
thanides are characterized by their exceptionally large 
electron orbital angular momentum, which lead to large 
anisotropic dispersion interactions between these atoms. 
Our theoretical estimate shows that in both Er and Dy 
collisions the ratio of anisotropic to isotropic dispersion 
interaction ACq/Cq is about 10%. This anisotropy leads 
to significant splittings among the 48 and 81 gerade 
short-range potentials that dissociate to the ground-state 
atomic limits of Er and Dy, respectively. We have shown 
that each potential has its own ro-vibrational structure, 
which by coriolis forces and the Zeeman interaction in¬ 
teracts with that of other potentials, creating a dense 
distribution of levels near the threshold and initiating 
chaos. In fact, we find a very large number of partial 
waves contributing to the creation of Fano-Feshbach res¬ 
onances. 

On the other hand, if we just consider the anisotropy 
from the magnetic dipole-dipole interaction alone, our 
coupled-channel calculations indicate that chaos in the 
level distribution does not appear. The strength of the 
dipole-dipole interaction is too small. In addition, we 
have shown that the NNS distributions for Dy and Er 
are very similar, as can be expected from their similar 
ACe/Cg ratio. The difference in their magnetic moment 
only plays a small role. This further confirms that chaos 
is due to the anisotropic dispersion interaction. 

The distribution of Feshbach resonances of ultra-cold 
ground-state alkali-metal, alkaline-earth, Yb, and Cr 
atoms, as experimental studies have shown, are not 
chaotic. This is because these atoms have a zero electron 
orbital angular momentum and, hence, only an isotropic 
dispersion interaction. Eventhough alkali-metal and Cr 
atoms have a non-zero magnetic moment of l/is and 6/i_B, 
respectively, these moments do not lead to chaos. We 
would expect that other magnetic lanthanides and ac¬ 
tinides with non-zero orbital angular momentum will ex¬ 
hibit chaotic properties in their collisions. In addition, 
collisions between mixed species, such as magnetic lan¬ 
thanides and alkali-metals, like K-fDy or Na-|-Er, might 
be susceptible to chaos. A first theoretical analysis for 
Li-|-Er [59], however, estimates a small 2% dispersion 


anisotropy and no chaos is predicted. 

Another interesting property of magnetic lanthanide 
gases is the extreme sensitivity of the atom-loss spec¬ 
tra and, in essence, three-body recombination to the 
temperature. This phenomenon was first observed in 
Ref. [22] for loss spectra of Dy. The number of Dy res¬ 
onances increases by 50% when the temperature is in¬ 
creased from 420 nK to 800 nK. Here, we observe a 25% 
increase in the Er resonance density when the temper¬ 
ature rises from 250 nK to 1400 nK. We have shown 
by a comparison of resonance profiles taken at several 
temperatures and predictions of a theoretical model of 
three-body recombination via the formation of a trimer, 
or more precisely of a shared pair bound state, that the 
origin of the temperature-dependent resonances lies in 
the “partial wave” of the three-atom entrance channel. 
Entrance channels with zero and non-zero total orbital 
angular momentum M lead to line shapes with a dif¬ 
ferent temperature behavior. Those with A/” = 0 or 
“s-wave” entrance channels have sharply decreasing re¬ 
combination rates with temperature, whereas those with 
A/" = 2 or “d-wave” entrance channels have an increasing 
recombination rate. Temperature sensitive resonances 
can only be explained by “d-wave” collisions. It is worth 
noting that for alkali-metal-atom collisions a number of 
entrance-channel p-wave resonances have been observed 
(See for example Ref. [50] for cesium). Analysis of the 
temperature dependent rate coefficient, however, was not 
performed. 
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